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We perform non-perturbative calculation of the tv° — > 77 transition form factor and the associated 
decay width using lattice QCD. The amplitude for two-photon final state, which is not an eigenstate 
of QCD, is extracted through an Euclidean time integral of the relevant three-point function. We 
utilize the all-to-all quark propagator technique to carry out this integral as well as to include 
the disconnected quark diagram contributions. The overlap fermion formulation is employed on the 
lattice to ensure exact chiral symmetry on the lattice. After examining various sources of systematic 
effects except for possible discretization effect, we obtain r„.o^ 77 = 7.83(31)(49) eV for the pion 
decay width, where the first error is statistical and the second is our estimate of the systematic 
error. 
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The neutral pion decay process provides a unique op- 
portunity to test a fundamental symmetry property of 
the gauge theory. A quantum effect due to a fermion 
loop violates the axial-current conservation, and gives 
the dominant contribution to the 7r° — > 77 decay 
rate. The prediction from this Adler-Bell-Jackiw (ABJ) 
anomaly [HH] (or the chiral anomaly) is rather precise be- 
cause higher-loop diagrams do not contribute in the limit 
of vanishing quark mass and external momenta [3] even 
under the presence of strong interaction. On the other 
hand, a recent experimental measurement of the neutral 
pion decay width has reached the accuracy of 2.8% [2] 
and is aiming to achieve 1.4% in the near future. At this 
level of accuracy the correction due to finite quark mass 
becomes relevant. Phenomenologically, an enhancement 
of the decay width of about 3-5% has been expected ac- 
cording to the sum rule and chiral perturbation theory 
(xPT) approaches [5 8J, in which the main effect comes 
from a mixing of ir° with 77 and 77' mesons. For a recent 
review, we refer the reader to Ref. [5]. 

In this letter we present a model-independent calcula- 
tion of the 7T° — > 77 amplitude using the lattice Quantum 
Chromodynamics (LQCD) including dynamical up, down 
and strange quarks. We use the overlap fermion formula- 
tion [10] , which preserves chiral symmetry at finite lattice 
spacings. In this formulation, the chiral anomaly appears 
through the Jacobian of chiral transformation; Atiyah- 
Singer's index theorem is reproduced as far as the back- 
ground gauge field is smooth enough llj. On the other 
hand, whether the chiral anomaly is correctly reproduced 
at practically used lattice spacings (~ 0.1 fm) is a non- 
perturbative problem, that we address in this work. 



Compared to previous attempts [T2J [T3] , a new tech- 
nique is applied to treat two-external-photon state, by 
utilizing the all-to-all quark propagator [TU [15] . In the 
limit of degenerate up and down quark masses we obtain 
the decay rate with a statistical error of 4% and a total 
error of 7% after examining possible systematic effects. 

The 7T° — > 77 decay rate at the leading order of QED 
can be expressed as 

7T { \ III 

r^o^— p^^mJ.O.O), (1) 

where a e is the fine structure constant, is the neutral 
pion mass and J r 7r o 77 (m 2 i s the form factor of 

the pion to two (virtual) photon transition with p\^ the 
photon momenta. In the chiral limit, the ABJ anomaly 
predicts 

= ^o 77 (0,0,0) = , (2) 

where Fq is the pion decay constant F n in the chiral limit. 
We define a normalized form factor as F(to 2 ,p 2 ,p 2 ) = 
(47T 2 i^)J r 7r o 77 (m 2 ,Pi,f>|)- In the Minkowski space-time 
J 7 w o 77 (to 2 .,p 2 ,P2) is defined through the matrix element 

= e^apPiP^F^o^ml^l^l) (3) 

where q is the 7r° momentum satisfying the on-shell con- 
dition q 2 = to 2 . The current j M = YlfQfQfliiQf ^ s 
the hadronic component of the electromagnetic vector 
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current and the sum runs over all relevant quark fla- 
vors: / = u,d,s. Qf denotes the electromagnetic charge 
of them: Q u = +2/3 and Qd :S = —1/3. The factor 



£ fj.ua/3P1P2 i s induced by the negative parity of 7r°. 

By an analytic continuation of ([3| from the Minkowski 
to Euclidean space-time [El [17] one may write 



M^ v {p u p 2 ) = lim 
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2£+ff 



dh e^-^C^ih^U) 



C^{hM^)= J d 3 5e-^^y' ( i 3 Fe^^O|T{i M (x,ti)j I/ (0,t 2 )7r (^^)}|O) , 
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where ti, £ 2 an d ^ are Euclidean time slices. 
/ d 3 z e lq ' z Tr°(z, t n ) is an interpolating operator for the 
neutral pion with the spatial momentum q. Its amplitude 
and energy in the ground state are denoted by (j)^^ and 
E v ff, respectively. The four-momentum of the first pho- 
ton pi = is chosen as input, while the momentum 
of the second photon is given as pi = (E Wi ^—lo, q—pi) by 
momentum conservation. Note that the analytical con- 
tinuation is valid only for p\ 2 < My, which sets the limit 
on the value of ui. (My stands for the invariant mass of 
the lowest energy state in the vector channel.) Since the 
two photons can not be on-shell simultaneously, we cal- 
culate the form factor at the off-shell photon momenta 
and then extrapolate to the on-shell limit. 

To calculate the matrix element (tt\j^jvir \Q) in 
Eq. @, we use 2+1-fiavor overlap fermion configura- 
tions generated by the JLQCD and TWQCD Collabo- 
rations [TH1 US] at a single lattice spacing a — 0.11 fm 
and two spatial lattice sizes L/a = 16 and 24. The time 
extent is T/a = 48. Although the main gauge ensem- 
bles have a fixed (global) topological charge Q = 0, the 
deviation from the 9- vacuum is understood as a finite vol- 
ume effect of 0(l/L 3 T) 20J . We check the significance 
of this effect by comparing the results with two different 
values Q = 0,1. We utilize the all-to-all propagator to 
calculate the correlation function C '^(ti, t2, i^) at any 
time slices of ti, i 2 and t^. The electromagnetic cur- 
rent jp_ — ^2fQfQf"tfiQf is implemented on the lattice as 
a local operator with a renormalization factor calculated 
nonperturbatively in |21] to match the lattice results with 
the continuum theory. The up and down quarks are de- 
generate in mass. We use the bare values am u .d — 0.015, 
0.025, 0.035 and 0.050, corresponding to the pion mass 
ranging from 290 to 540 MeV. Our final results are 
obtained by an extrapolation of the data to the the phys- 
ical pion mass m T p hy The strange quark mass is fixed 
at am s = 0.080, which is very close to the estimated 
physical value. 

From the large tip — behavior of C '/^(ii, i 2 , i^), it 
is possible to extract the 7r°-ground state. We define the 
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FIG. 1: The amplitude A n (r) as a function of r for momen- 
tum setups #1 (left) and #2 (right). The black (red) curves 
indicate the lattice (VMD) amplitudes. 



amplitude A v as 

K{r)= lim C^ituh^/e-^*-^ 



(5) 



with r = ti — ti and t = rmn{ii,i 2 }, and obtain 
^■1[iv{pitP2) by performing an integral 



2-Krr a 



dre ur A n ( T ) + / dr e^- E ^ )T A„(T) 
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We use two momentum setups p\ = ^(0,0,0), q — 
f (0,0,1) (setup #1) and pi = ^(0,0,1), q = 
^(0,0,0) (setup #2). The resulting amplitudes A n (r) 
for these setups are shown in Fig. [I] In order to qual- 
itatively understand the r-dependence of the A 1t (t) 1 
we consider the vector-meson-dominance (VMD) model 
^°K,pf,pl) = c v G v (pI)G v (pI), with Gy(p 2 ) = 
My /(My — p 2 ) the vector meson propagator and cy a 
constant. The amplitude A% mu (t), reconstructed from 
this model, is plotted by red curves in Fig. [l] (The de- 
tailed expression for A^ md (t) will be given in a later 
publication [22].) We find that the VMD model describes 
the lattice data already at \r\/a = 7, and we can safely 
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FIG. 2: Left: Contour of {p\,v\) rescaled by 1/M V for our 
momentum setups. Right: F 7r o 77 (m^.,p 1 ,p2) as a function of 
p\/My. Lattice data are obtained on a 24 3 x 48 lattice at 
am u ,d = 0.015. 



evaluate the contribution beyond \r\/a = 13, where the 
lattice data are truncated due to the finite time extent 
T. At small |r| the VMD model fails to match the lattice 
data. This is because no information of the vector-meson 
excited states are contained in J^oH^ \m\,p\,p\^). Given 
the dominant role played by the lowest vector meson, we 
take its form as a basis to analyze the functional form of 
T 7r o 11 {m'^,p\,p\). Namely, we use an expression 



•^77 ("4 > Pi > Pa) = c vGv(pl)G v {pl) 

+E c - ((pirovipi) + ( P i) m G V (p 2 2 )) 



(7) 



which includes possible contributions from excited states 
as a polynomial of p\ 2 . In the chiral and large volume 
limit for which the two-pion threshold opens, the VMD 
model would no longer give an adequate description of 
the data. Such effect could become significant for pre- 
cision better than few percent, which is not within the 
scope of this work. 

By varying u, we obtain M fll/ (pi,p2) in a certain range 
of p\ and p\. As shown in the left panel of Fig. [2j a pair 
Oi.Pl) = -Pi> ( E *,q ~ w ) 2 - (<f-.Pi) 2 ) forms a con- 
tinuous contour on the (pi,pf) plane for p\ 2 < My/2. 
Evaluating F 7t a 1 ^ l {m'^ 1 p\ 1 p\) along this contour, we ob- 
tain the data plotted in the right panel of Fig. 2] We per- 
form the combined fit of these data to Eq. ((7) with four 
free parameters: cv, Cq, and 00,1 = c^o, truncating 
the higher-order terms which turned out to be negligibly 
small. The fitting curves are shown in the right panel of 
Fig. [2] As expected, the single formula describes the 
data with different momentum setups. Combining the 
resulting fit parameters, we obtain the normalized form 
factors F(m%,0, 0), which are plotted in the uppermost 
panel of Fig. [3] 



In the following we analyze the details of systematic 
effects. When calculating the integral in Eq. ([6]), we use 
the summation instead of the integration. This causes 
a discretization effect, which vanishes in the continuum 
limit. Putting A^ md (t) into Eq. we find that the 
fractional difference between M PiV (j>\,p2) from the sum- 
mation and the integration is less than 5 x 10~ 4 . With 
the lattice data that include the excited state contribu- 
tions, we could expect a larger error, ~ 1 x 10~ 3 , which 
is estimated from a difference between VMD and lattice 
data in Fig. [T] We can therefore safely neglect this source 
of error as it is well below 1%. 

We use two lattice volumes and two topological- 
charge sectors to check finite-size (FS) effects. Follow- 
ing Ref. [2U] we analyze the fixed-topology (FT) effect 
and find it suppressed due to the kinematical structure 
of £^vapPiP2- By comparing the lattice results at dif- 
ferent topological-charge sectors, we do not observe sta- 
tistically significant FT effects. The leading FS effect in 
G^ v {t\, t2, t v ) is the conventional one and known to be- 
have as e~ mxL 23 . To reduce the contamination due to 
this effect, we therefore use the data with m^L > 4 to 
perform the chiral extrapolation. (Namely, we exclude 
the L/a = 16 data points at the lowest two pion masses.) 

%PT shows that up to next-to-leading order (NLO) 
the m„.-dependence of -F(m^.,0,0) involves no chiral log- 
arithm [25]. We therefore simply fit F{m\, 0, 0) by a 
linear function in m^., and obtain F(0,0,0) = 1.016(20) 
and F(ml phy , 0,0) = 1.011(19). To check the higher- 
order correction, we also perform a quadratic fit under 
the constraint from the ABJ anomaly: F(0,0, 0) = 1. 
We do not find any statistically significant difference due 
to the higher-order term. The linear (quadratic) fit is 
shown by the solid (dashed) line in the uppermost panel 
of Fig. [3] 

Next we consider the data with m^L < 4, which tend 
to suffer from the FS effect. As shown in Fig. [3j at 
uift « 290 MeV we find that F(ml,0, 0) calculated at 
L/a = 16 lattice is 27% less than the one at L/a = 24. 
Although large, such FS effect is understandable. By in- 
serting the ground state into (j^ji/7r°), we can approx- 
imate this three-point correlation function with three 
hadronic matrix elements: 

(j^n ) -> (n\j,\V,e)(V,e\Un )(n \7T°\n) . 

The first matrix element is related to the electromag- 
netic coupling gv as (n\j^\V,e) = MygyE^, the sec- 
ond is proportional to the Viry coupling fly™ and the 
third is related to F„ by the PCAC relation. In our cal- 
culation we do not observe significant FS effect in My 
but find 8%, 7% and 9% shifts in gv, gv-n-y and F n , re- 
spectively, from L/a = 16 to 24, as shown in Fig. [3] 
These FS effects may accumulate in the three-point func- 
tion. We estimate the FS corrections Rq = 0(oo)/0(L) 
with O = gv, gv-n-y and F n . R gv and R gv are evalu- 
ated by adding a correction term, e _m,ri , into the lin- 
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FIG. 3: F(m^, 0,0), gv, gvn-i, F v and FS corrected 
F(m^, 0,0) as a function of m% from top to bottom pan- 
els. In each panel, data with (L/a,Q) = (16,0), (24,0) and 
(16, 1) are plotted by the blue, red and green symbols, respec- 
tively. The yellow symbols indicate the Particle Data Group 
(PDG) [26 or PrimEx [3] experimental values for reference. 
The solid (dashed) curves show the result of the fit to the 
linear (quadratic) function. The dataset used in the fit is 
explained in the text. 
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FIG. 4: A comparison of the full contribution to F^m^O, 0) 
and the connected-only piece. The upper (lower) panel shows 
the result without (with) FS correction. 



ear fit form in the chiral extrapolation of each quantity. 
With such corrections taken into account we confirm that 
their chiral limit is consistent with experimental data. 
R Fir is calculated to NNLO by using %PT 27 j. As- 
suming that i?F(m 2 , o,o) = RgyRgv-n-yRF* we ma y correct 
F(m^,0, 0) by a factor of RF(m 2 , o.o)- As shown in the 
lowest panel of Fig. [3] with FS correction F(m^., 0, 0) at 
L/a = 16 agrees with those at L/a — 24. Using the cor- 
rected data to perform a linear extrapolation, we obtain 
F(0,0,0) = 1.045(35) and F(ml phy ,0,0) = 1.041(32). 
The difference between the results from the two methods 
is considered as a systematic error. 

So far, our results are obtained neglecting the effect 



of disconnected diagrams that may appear because the 
electromagnetic current j M contains flavor-singlet con- 
tribution. Calculation of the disconnected diagram is 
computationally demanding and statistically noisy. We 
solve these problems by the use of the all-to-all prop- 
agator. The full data, including both the connected 
and disconnected contributions, are plotted in the up- 
per (lower) panel of Fig. [1] for the case without (with) 
the FS correction. We find that, although not signifi- 
cant, there is a shift from the connected data to the full 
ones. Since the accuracy of our calculation reaches a few- 
percent level, the disconnected effect is relevant. Using 
the full data we repeat the analysis. The linear fit of 
F(ml, 0,0) with m^L > 4 yields F(0,0,0) = 1.009(22) 
and F(m^ phy! 0,0) = 1.005(20). The fit with FS cor- 
rected F(ml, 0,0) produces F(0,0,0) = 1.007(36) and 
F(ml iphy , 0, 0) = 1.006(33). Including the disconnected 
contributions, the normalized form factor in the chiral 
limit and at the physical pion mass shifts by 1-4%. This 
is comparable to the statistical error. 

Using the full data we quote our results for F{m^, 0, 0) 
and r 7r o 77 in the isospin symmetric limit as 

F(0,0,0) = 1.009(22)(29) , 
FK phy ,0,0) = 1.005(20)(30) , 

r> 77 = 7.83(31)(49) eV . (8) 

where the systematic errors originate from the difference 
of the results by using two methods of treating the FS ef- 
fect. (The difference appearing in the full data is small. 
To be conservative, we use the connected data to esti- 
mate such systematic error.) Our results reproduce the 
predication of the ABJ anomaly F(0, 0, 0) — 1 and agree 
with the PrimEx measurement r 7r o 77 = 7.82(22) eV [4]. 
For future improvements, isospin breaking effects due to 
the light quark mass difference need to be included. 
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